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Abstract 

Within the next few years experiments at RHIC and the LHC wih seek to create in the laboratory 
a quark-gluon plasma, the phase of matter in which the Universe was initially created. It is believed 
that the plasma will survive long enough to reach thermal equilibrium. I give an introduction to the 
formalism of thermal field theory, the combination of statistical mechanics and quantum field theory 
needed to describe the plasma in thermal equilibrium, in a way that tries to keep close to the physics 
it describes. 



Introduction 

Thermal field theory is a combination of quantum field theory and statistical mechanics. This means 
that it is both difficult and interesting. The reason that we study it is that we want to describe the 
quark-gluon plasma, the phase that matter is believed to take above some critical temperature T^. 
Lattice calculations suggesft^l that Tc is about 100 MeV, or lO^^K. In the plasma phase the quarks 
and gluons are deconfined; they can move rather freely through the whole plasma. This is the phase 
in which the universe was created at the big bang, and before the end of the century experiments at 
the new collider RHIC will try to re-create it in the laboratory, by making gold nuclei collide together 
head-on and dump their kinetic energy into a small volume. Similar experiments, at much higher 
energy, are planned later for the LHC at CERN. 

There is an obvious question: if a plasma is indeed produced, how will we know it? As yet there is no 
simple answer. There are estimates t^l , necessarily based on very crude non-equilibrium theory, that 
suggest that the plasma will survive for a time long enough that it reaches thermal equilibrium before 
it eventually decays back into ordinary matter. So far, it is only equilibrium thermal field theory that 
is well formulated, and my lectures concentrate on this. For more information, there is a book that 
was published last yeart"^] and is already the standard text. As I want my description to stay as close 
as possible to physics I will develop the theory using operators rather than path integrals, and mostly 
I will use the so-called real-time formalism. 

Because in relativistic theory particles are continually being created and destroyed, it is appropriate 
to use the grand partition function 



Z = j;(z|e-^(^-^^)|i) (1) 

i 

Here (3 is the inverse temperature, (3 = l/ksT, and usually we use units in which Boltzmann's constant 
= 1. The system's Hamiltonian is H and is some conserved quantum number, such as baryon 
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number, with the corresponding chemical potential. There may be several conserved quantum 
mumbers, in which case is replaced with y^/jgiVQ,. 

a 

The states \i) are a complete orthonormal set of physical states of the system. In scalar field theory 
all states are physical and so 



Z = tr e-^^^-"^) (2) 

which is invariant under changes in the choice of orthonormal basis of states. In the case of gauge 
theories there are unphysical states, for example longitudinally-polarised photons or gluons, which 
must be excluded from the summation in (1). So then 



Z = tr P e-P^^-"^^ (3) 

where P is a projection operator onto physical states. The presence of P can make things more 
complicated, and so to begin with I will consider scalar filed theory, where it is not needed. 

All the macroscopic properties of the system in thermal equilibrium may be calculated from Z. In 
particular, for a system that is so large that its surface energy is negligible compared with its volume 
energy, the equation of state is 



PV = T log Z (4) 
Also, the "thermal average" of an observable corresponding to an operator Q is 



<Q>= Z-^ tr Qe-^^^-i"^^ (5) 

Notice that, throughout, all operators are familiar zero-temperature ones. The temperature enters only 
in the exponential, which characterises the particular ensemble of states used to calculate expectation 
values of the operators. 



Noninteracting scalar bosons 

For many systems of bosons there is no conserved quantum number N ; for example, in the case of a 
heat bath of photons there is no constraint on their total number. Then in the scalar-field-theory case 
Z is just tr e~^^ . 

In the absence of interactions, the energies of the separate particles arc good quantum numbers. 
To begin with, quantise the system in a finite volume V, so that the single-boson energies er are 
discrete. The states |i) of the system are labelled by the single-particle occupation numbers rir, and 
the eigenvalues of the noninteracting Hamiltonian Hg are 

niei + n2e2 + wsCs + • ■ ■ 
So the noninteracting grand partition function is 

_ ^ g/3(niei+n2e2+...) 

{rir-} 

r rir r 
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and so 



log^o = -i;iog(l-e-^^0 (66) 



In the continuum limit 

(27r) 

and so the nonintcracting equation of state is 



F^flogZo^-ry^loga-e-*) (8) 

where kg = Vk^ + fn^- (If the bosons have non zero spin, there is an additional factor gg corresponding 
to the spin degeneracy of each single-boson state.) 

In the continuum limit we usually work with fields 
with 

[a(k),at(k')] = {2Trf2koS^^^k-k.') (10a) 
In the discrete case, we usually define 

[ar,al] = 5rs (106) 

If we sum this over r, the result is 1. But if we apply V J d^k/(2'K)^ to (10a), the result is rather 
2k^V . That is (10a) and (10b) have definitions of the operators a differing by a factor '^2koV . We 
correct for this by defining the field in the discrete case to be 

(t){x) = ^ -^=L=a^e-*^'"*e^'''^-'^ + h.c. (96) 

We can now calculate the thermal average 

{Tct>{x)m)o = ^o"' l;(^|e-^^''^0(x)<^(O)|^) (11) 

i 

Observe first that 

Tct,{x)m = {o\Tct>{x)m\o)+ ■■ Hx)m ■ (12) 

where : : denotes the usual normal product. The first term contributes to (11) just the usual zero- 
temperature Feynman propagator. To evaluate the contribution from the second, use the discrete case 
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(9b) and so obtain double sums of terms i . When we take the necessary expectation 

r,s 

values, only the first survives, and then only for r = s. We may replace aja^ with n^, and use 

^ ^ ' Y,ne--^^ = f{e) (13a) 



1 - e-P^ 



where / is the Bose distribution 



f{e) = (136) 



So we find that 



(: mm ■.)o = ^Yl ^/(er)e*^'^*-^'''^-^ + h.c. (14) 



V ^ 2er 

r 



Going to the continuum limit, we have 



mx)m)o = / e-^-i?o (fc) 



(15) 



where n{k'^) = /(|A;°|). The second term is the contribution from the heat bath; it contains the 
(5-function because so far the heat-bath particles do not interact and so they are on shell. 

Perturbation theory 

Suppose now that we introduce an interaction and again calculate {T(/){x)(/){0)) . Now 4>{x) is the 
interacting Heisenberg-picture field: it is the familiar operator of zero-temperature field theory. We 
can develop a perturbation theory along the same lines as at zero temperature, by introducing the 
interaction picture that coincides with the Heisenberg picture at some time to: 

</.j(i,x) = A(i)</>(t,x)A-i(t) 

^ ^i{t-to)Hoi ^-i{t-to)H V ) 

where Hqi is the free-field Hamiltonian in the interaction picture. As usual 

(17) 



Texp (^-i J^' dtHf^it)^ 



We need U{ti,t2) for complex ti and t2, so that we integrate t along some contour running from t2 to 
ti in the complex plane and generalise the time ordering T to an ordering T^. along C: The operator 
whose argument is nearest to ti along the contour comes first. 

Now, from the definition of U, 

^-0H ^ ^-pHoi u{to-if3,to) 
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so that 



Z-hi e-^-^(^(x)(^(0) = Z-hr e-^^°'U{to - iP, x°) (/)j(a;)[/(ar°, 0) (f)i{0) U{0, to) 

= ZoZ-^U{to - iP, x°) ct)i{x) C/(xO, 0)0/(0) U{0, to))^ ^"^^^ 
where I have used U{ti,t2)U{t2,t3) = U{ti,t3). 
Compare with what one has at zero temperature: 

(O|0(x)(^(O)|O) = (O|C/(oo,x°)07(x)[/(x°,O)(^7(O)C/(O,-oo)|O) 

In (18) we have a non-interacting thermal average instead of a vacuum expectation value, {to — i/3) 
instead of oo, and to instead of — oo. In fact there is a close similarity between thermal perturbation 
theory and the usual zero-temperature Feynman perturbation theory. The main difference is that 
instead of the internal lines in the ordinary Feynman graphs representing vacuum expectation values, 
in the thermal graphs they are non-interacting thermal averages. In order to derive this, one needs to 
establish Wick's theorem. It is a remarkable fact that indeed, for example, 

{n(l)i{xi)(t)i{x2)(l)i{x3)(l)i{x4))o = ^{T(t>i(t>i)o {T(t>i(t>i)o (19) 

where the sum is over the possible pairings of the fields. For almost every ensemble other than one in 
thermal equilibrium there would be correction terms to (19). 

One needs to choose a value for to- A common choice is to = 0, with the contour C for the t integrations 
running from to —i^ along the imaginary axis. This is the imaginary-time formalism. Alternatively, 
to — ^ — oo, which with suitable contour choice gives the real-time formalism. 

Real-time formalism 

One is interested in equilibrium properties of the plasma at finite times. Presumably these are inde- 
pendent of how it reached thermal equilibrium. So, as in familiar scattering theory, we are free to 
imagine that the interaction slowly switches off as we go into the remote past, and then when we take 
to — oo the interaction-picture fields become the usual noninteracting in fields. The corresponding 
in states are direct products of non-interacting single-particle states. We need to choose how the 
contour C runs from — oo to (— oo — i/3), and the choice that keeps the formalism in the most direct 
contact with the physics is the so-called Keldysh one: along the real axis from — oo to oo, back to 
— oo, then straight down to (— oo — i/?). 

V 



For most applications (though not alll^] ) it turns out that the vertical part of the contour does not 
contribute. Then we may write the propagator that corresponds to a line of a thermal graph as a 2 x 2 
matrix: 



D(xi,a;2) 



(T0in(xi)(/)i„(a;2))o 

{(t>iAxi)4>ir,{x2))o 



{x2)4>\n{xi))o 

(T0in(a;i)0i„(x2))o 



(20) 



When both xi and X2 are on the — oo to oo part of C, the ordering is ordinary time ordering 
T; this corresponds to the element Dn of D. When both are on the oo to — oo part of C, is 
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anti-time-ordering T; this corresponds to -D22- The off-diagonal elements correspond to xi being on 
one part of C and X2 on the other. There is translation invariance: the elements of D depend only on 
the difference between xi and X2- 

I have already shown how to calculate Dn; the result is given in (15). The other elements of D may 
be calculated in the same way, and its Fourier transform is 



k'^—m^+ie 
U2 



It may also be written in the form 



27r<5-(A;2-m2) 



(21a) 



M D M 



(216) 



with 



D =i 



k'^—m'-'+ie 





-1 



k^ — —ie 





e■2^ 

g 2 



62 



/3fc" 



\p\k°\ 



e 2 



(21c) 



For the case of a fermion field, there is a rather similar matrix propagator, but with the Fermi-Dirac 
distribution replacing the Bose distribution. 



Matrix structure 

The elements of the matrix propagator (20) are not independent. For example, 

D2i{x) = (</.i„(x)(/)i„(0))o = Z^hv e-^^°''>(/)in(2:)</)i„(0) 



[X e 



= Z^\t e-^^°-(/>in(0)e 
= Zo"Hr e-^^"-(/)(0) (/)(x° - i/3,x) 
= L»i2(x°-i/3,x) 



/3//oin 



(22a) 



Here, I have used a general property of traces, that iT:{AB) = tr(i?^), and the fact that Hqi^ is the 
time-translation operator for the noninteracting field (?!)in. The Fourier transform of (22a) is 



D2i{k) 



(226) 



Also, from their definitions (20), one can see that Dn and D22 may be expressed in terms of D12 and 
D21. It is this, together with (22b), which is responsible for the matrix structure (21). 

Define a dressed thermal propagator matrix D'(xi,X2) analogous to D(xi,a;2) in (20), but with the 
interacting Heisenberg field instead of 0in. For example, D'^2{^'i.^^2) = {(f>{x2)(t>{xi))- Then, because 
H is the time-translation operator for 0, we can again derive 



(22c) 



and so deduce that D' has the structure 
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(23) 



Define the thermal self-energy matrix 11 by 



(24a) 



Then 11 has the structure 



,n = M-r^n'(fc,r) 



\iU'ik,T)]* 



(246) 



If we then solve (24a) for D', we find 







fc2-m2-n* 



M 



(25) 



So it is natural to interpret Re ft as a temperature-dependent shift to the mass m^. ft also has an 
imaginary part, so the propagation of the field through the heat bath decays with time. 

In scalar field theory, 



n = 





To calculate the contribution to II12 from the second term, for example, one needs 




2 + 




2 + 1 




2 + n 



where each line r 
opposite in sign. 



1 represents Dij[k) and each vertex 2 is the same as the normal vertex 1, but 



Imaginary-time formalism 

The real-time formalism stays close to the physics, but has the calculational complication that the 
propagator is a matrix. In the imaginary-time formalism there is not this complication, though except 
for a few simple cases there is the need to perform an anlytic continuation from imaginary to real 
time at the end of the calculation. In the imaginary-time formalism the t integration runs along the 
imaginary axis, so ordinary time-ordering is replaced with ordering in imaginary time: 



D{xi,X2) = ^(— Im t)D2i{xi,X2) + 0{lm t)Di2{xi,X2) 



(26) 



where t = — X2- Because both and are integrated from to i(3, we need D{xi,X2) for values 
of Im t in the range —f3 to +p. In this finite interval it has a Fourier-series expansion. 
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D{t,^) = ^ E ^n(x)e-'^* (27) 

n=— oo 

where lj„ = nvr//?. However, the relation (22a) impHes that D{t,x) = D{t + so that only 

even values of n contribute to the sum. (In the case of fermions, the anticommutativity of the fields 
results in a minus sign appearing in the corresponding relation (22a), and so then only odd values of 
n contribute.) 

If we apply a 3-dimensional Fourier transformation to (27) and invert the Fourier summation over n, 
we find 

L>„(k)= / dte-«'"*L>2i(t,k) (28) 

which turns out to be just the ordinary Feynman propagator with = itOn- So the Fcynnian rules 
are just like the zero-temperature ones, except that the energy-conserving (5-function at each vertex 
is replaced with a Konecker delta which imposes conservation of the discrete energy, and round each 
loop of a thermal graph 

d-^k f d^k , , 



(2;r)4 (27r 



n 



Gauge theories 

For gauge theories there is the complication that the grand partition function has to include the projec- 
tion operator P onto physical states: see (3). There are two formalisms for the resulting perturbation 
theory!^] : 

A Only the two physical degrees of freedom of the gauge field (the transverse polarisations) acquire 
the additional thermal propagator; the other components of the gauge field, and the ghosts, 
remain frozen at zero temperature. (This is for the bare propagators; self-energy insertions in the 
unphysical bare propagators do depend on the temperature.) 

B All components of the gauge field, and the ghosts, become heated to temperature T . 

In the zero-temperature field theory, the ghosts are introduced in order to cancel unwanted con- 
tributions from the unphysical components of the gauge field, and the two formalisms lead to the 
same answers for calculations of physical quantities for that reason. Often, using formalism A makes 
calculations simpler. It also makes them stay closer to the physics. 

Photon or dilepton emission from a plasma 

As an application, consider the emission of a real or virtual photon of momentum q from a quark-gluon 
plasma. This is supposed to be an important diagnostic test of whether a plasma has been created in 
an experiment and reached thermal equilibrium, and a way to measure its temperature. 

Before a photon is emitted, the plasma is described by the density matrix 

p = .Z-i^|iin)(iin|e-^-^ (30) 

i 

The emission probability is calculated from squared matrix elements of the Heisenberg-picture elec- 
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tromagnetic current: 

W^iq) = Z-^Y^ J d'^xe-"i-''{fout\J^'{x) |un)(wn|e-^^ j J'^(0)|/out) (31a) 

where I have introduced also a complete set of out states for the plasma. These satisfy the completeness 
relation 



^ I/out) (/out| 

/ 



and so 



W^'^'ig) = Z-^Y, J (^^a;e'«-^(iin|e-^-f^J^(0)J^(x)|iin) (316) 



If we introduce a matrix G^'^(q') analagous to D', but with thermal averages of products of electro- 
magnetic currents instead of fields, W'^'^{q) is just G^'j (—<?)• So we draw thermal graphs where current 
—q enters at a 1 vertex and leaves at a 2 vertex, and distribute the labels 1 and 2 in all possible ways 
on the other vertices. For example. 




The emission rate is calculated from a matrix element times its complex conjugate; the vertices labelled 
1 correspond to a contribution to the matrix element and those labelled 2 to its complex conjugate. 
The sets of 1-vertices and of 2-vertices are joined by 12 lines which, according to (21a) are on shell 
and represent particles in the heat bath. The thermal graphs sum together many physical processes. 
Consider the first graph, for example. Its right-hand part represents the contributions 




to the amplitude. 

In fact, energy-momentum conservation allows only the first one to be non-zero. I have not drawn in 
the other heat-bath particles, but remember that they are there as spectators. The other part of the 
thermal graph contains only 11 lines. If I use only the zero-temperature part of Dn in each, I obtain 
the amplitude 



(plus other terms which again vanish for kinematic reasons) and so part of the thermal graph represents 
the interference between this and . If instead I use the thermal part n{k^) 2Tr6(k'^) of the 11 gluon 
propagator, I obtain the amplitudes 
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In each case, the incoming and outgoing gluon Hnes must have the same momentum k, so that these 
amphtudes again interefere with ^^^^^ , but now with the gluon k being one of the spectator particles 
in the heat bath: 




k 



Similarly, I can identify physical processes that involve the thermal parts of the ii quark propagators. 

Even a simple-looking thermal graph corresponds to a large number of physical processes t^l, each of 
which can be rather complicated. An example is 




2 2 2 

for which just one of the physical processes is the interference between 




A disconnected graph occurs because the thermal graph has an "island" of a pair of 1-vertices entirely 
surrounded by 2-vertices 

Infrared divergences 

The infrared divergences of zero-temperature become much worse at finite temperature: the Bose 
distribution diverges at zero energy and causes the usual logarithmic divergences to become power 
divergences. We know that the infrared divergences must cancel if the theory is to make sense, and 
in practice they always do, but there is no general theory to show this. To some extent, the situation 
can be rescued by including thermal self-energy insertions in the propagators, so that they acquire a 
mass proportional to the temperature. But, in the case of photons or gluons, not all the degrees of 
freedom have a mass, according to perturbation theory. If the so-called magnetic mass is non zero, it 
is nonperturbative. 
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Consider, for example, the effect on the decay rate vr — > e^e of the microwave background, which 
is a heat bath consisting only of photonst^L Let F be the decay rate in vacuum. The heat bath will 
change it partly because it gives the electrons an additional temperature-dependent mass Sm^l oc e^T^. 
This causes a change ST = 5ml dV/dml, which is associated with thermal graphs of the form 




There is also the thermal graph 




One finds that 



AT 



where 



6ml ^ m^^aEM 
7r3Q 



dmi 



+ - 



(32) 



In the integral, the first c^'-'*-' -function corresponds to the contribution from the internal vertices in the 
thermal graphs being both i or both 2 and the second to 12 and 21. Each term separately is infrared 
divergent, like J" dk/k"^, but the divergences cancel. 

In fact there is more cancellation than just that of the infrared divergences. For T <C Q one can 
expand (32) in powers of /Q^. One finds that the first term, of order OiEMT"^ /Q"^ ■, exactly cancels 
the electron- mass-shift contribution 5T, so that the net change in the decay rate is of order uemT^. 
The lesson is that, in thermal field theory, it is of importance to calculate all terms; gauge theories at 
finite temperature are rife with cancellations. 

Another example where there is an apparent infrared problem is that of the calculation of the equation 
of state for the quark-gluon plasma. For instance, in the purely gluonic case thermal graphs of the 
form 



- " r 



become more and more divergent as more vertical lines are added. To see this, use the imaginary time 
formalism. The term in the multiple sum over the energies for which all the energies vanish looks just 
like an integral that would occur in zero-temperature 3-dimensional QCD, and simple power counting 
at each k = reveals the problem. However, it seems that the problem goes away if one sums over all 
thermal graphs. Give each gluon a mass m. Then, apart from possible difficulties with doing this in 
a gauge theory, in the real-time formalism one can derive the formulate] 



(33) 
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where 11 is the self energy of the gluon defined as in (24b) and e{qo) = ±1 according to whether go 
is positive or negative. As any divergence of ft now appears in the denominator, the summation has 
made it harmless. One has to worry about the denominator possibly vanishing when q = and m = 0, 
but this will be rendered harmless by the appearing in d'^q = q^dqdil. 

Linear response theory 

Suppose that the thermal equilibrium of a plasma is disturbed by the switching on at i = of an 
external electrostatic potential A'^^^ {x) . Then the system's Hamiltonian acquires an extra term 

H'{t) = 9{t) I d'x Mx)A'',^{x) (34) 

where Jq is the charge density. This will cause Jq to change. As it is a Heisenberg-picture operator, 
its equation of motion is 

^=i[H + H',Jo] (35a) 

where H is the original Hamiltonian. When we take the thermal average of this equation, the con- 
tribution from H will disappear because originally there was thermal equilibrium. So the integrated 
change in {Jq{x)) at very large time is 

S{Mx)) = J d^x' Gr{x - x')Al^,{x') (356) 

where 

Gr{x - x') = e{t - t'){[Mx), Mx')]) (35c) 

Taking the Fourier transform, 

5{Mk)) = GR{k)Al^,{k) (35d) 
We may express the retarded Green's function Gr in terms of elements of the matrix Green's function 

G« = i(G;?-G°° + G°?-G°°) (36) 

Each of the terms here may be calculated from pertubation theory. However, it is simpler to express 
Gr in terms of the function &^{k) that appears in the diagonal matrix associated with (see 
(23)): 

Gr = n(fe°)e3^l'="l{(G°° + G°°*) sinh + (G°° - G°°*) sinh (37a) 

On the other hand 

Gn = n(feO)e^/3l'=°l{(GOO + (7°°*) sinh + (G°° - G°°*) cosh i/^fe"} 

So it is sufficient to calculate G^ and change its imaginary part by multiplying it by tanh ^pk^\ 
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